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Abstract 



In this paper we quantitatively discuss the impact of the current 
uncertainties in the even zonal harmonic coefficients of the Newto- 
nian part of the terrestrial gravitational potential on the measurement 
of the general relativistic Lense-Thirring effect using a suitable linear 
combination of the nodes f2 of the laser-ranged LAGEOS and LAGEOS 
II satellites. The 1-sigma systematic error due to the mismodelling in 
the J £ coefficients ranges from - 4% for the EIGEN-GRACE02S grav- 
ity field model to ~ 9% for the GGM02S model. Another important 
source of systematic error of gravitational origin is represented by the 
secular variations Jg of the even zonal harmonics. While the relativistic 
and Ji signals are linear in time, the shift due to Ji is quadratic. We 
quantitatively assess their impact on the measurement of the Lense- 
Thirring effect with numerical simulations obtaining a 10—20% 1-sigma 
total error over 11 years for EIGEN-GRACE02S. Ciufolini and Pavlis, 
in a test performed in 2004, claim a total error of 5% at 1-sigma level. 

Keywords: Lense-Thirring effect; LAGEOS satellites; GRACE Earth 
gravity field models; even zonal harmonics and their secular variations 



1 Introduction 

1.1 The Lense-Thirring effect 

The post-Newtonian Lense-Thirring effect (Lense and Thirring 1918, Soffel 
1989, Lammezahl and Neugebauer 2001) is one of the few predictions of 
the Einsteinian General Relativity Theory (GRT) for which a direct and 
undisputable test is not yet available. 

According to Einstein, the action of the gravitational potential U of 
a given distribution of mass-energy is described by the metric coefficients 
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9iivi fi, v = 0, 1, 2, 3 of the space-time metric tensor. They are determined, 
in principle, by solving the fully non-linear field equations of GRT for such 
a mass-energy content. These equations can be linearized in the weak-field 
(U/c 2 « 1, where c is the speed of light in vacuum) and slow-motion 
(v/c << 1) approximation (Mashhoon 2001; Ruggiero and Tartaglia 2002), 
valid throughout the Solar System, and look like the equations of the linear 
Maxwellian electromagnetism. Among other things, a noncentral, Lorentz- 
like force 

F LT = -2m (^) x B g (1) 

acts on a moving test particle of mass m. It is induced by the post- 
Newtonian component B g of the gravitational field in which the particle 
moves with velocity v. B g is related to the mass currents of the source of the 
gravitational field and comes from the off-diagonal components goi,i = 1, 2, 3 
of the metric tensor. Thanks to such an analogy, the ensemble of the grav- 
itational effects induced by mass displacements is also named gravitomag- 
netism. Far from a central rotating body of proper angular momentum L 
the gravitomagnetic field is 

G[3r(r.L)-r 2 L] 
B 9 = — 5 • ( 2 ) 

One of the consequences of Eq. (JJ) and Eq. @ is a gravitational spin- 
orbit coupling. Indeed, if we consider the orbital motion of a particle in the 
gravitational field of a central spinning mass, it turns out that the orbital 
angular momentum I of the particle undergoes the Lense-Thirring preces- 
sion, so that the longitude of the ascending node f2 and the argument of 
pericentre oj of the orbit of the test particle are affected by tiny secular 
rates ^lt, <^lt (Lense and Thirring 1918, Cugusi and Proverbio 1978, Soffel 
1989, Ashby and Allison 1993, Iorio 2001a) 

■ 2GL . 6GL cos i . 

"LT = r, ^LT = r> (3) 

c 2 a 3 (l — e 2 )a c 2 a 3 (l — e 2 )2 

where a, e and i are the semimajor axis, the eccentricity and the inclina- 
tion, respectively, of the orbit and G is the Newtonian gravitational con- 
stant. Note that in their original paper Lense and Thirring (1918) used the 
longitude of the pericentre w. 

In April 2004 the GP-B mission (Everitt et al. 2001) has been launched. 
One of its goals is the measurement of another effect induced by the ter- 
restrial angular momentum, i.e. the gravitomagnetic precession of the spins 
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(Schiff 1960) of four superconducting gyroscopes carried onboard with an 
expected accuracy of 1% or better. 

However, according to Nordtvedt (2003), the multi-decade analysis of 
the Moon's orbit by means of the Lunar Laser Ranging (LLR) technique 
yields a comprehensive test of the various parts of order 0(c~ 2 ) of the 
post-Newtonian equation of motion. The existence of the Lense-Thirring 
signature as predicted by GRT would, then, be indirectly inferred from the 
high accuracy of the lunar orbital reconstruction. Also the radial motion 
of the LAGEOS satellite would yield another indirect confirmation of the 
existence of the Lense-Thirring effect (Nordtvedt 1988). 

1.2 Aim of the paper 

In the case of the Earth, the precessions of Eq. Q are very tiny: for 
the laser-ranged LAGEOS satellites they amount to a few tens of milliarc- 
seconds per year (mas yr _1 in the following). Extracting such a minute 
signal from the background of the much larger competing classical effects of 
gravitational (even zonal harmonics and their secular variations, tides) and 
non-gravitational (direct solar radiation pressure, Earth albedo, direct Earth 
infrared radiation, thermal forces like the solar Yarkovsky-Schach effect and 
the terrestrial Yarkovsky-Rubincam effect) origin is a very challenging task 
that requires detailed knowledge also from many branches like geodesy geo- 
dynamics and celestial mechanics other than relativistic physics. In turn, 
the extremely high accuracy required for such relativistic tests in the ter- 
restrial space environment can help in further increasing our knowledge of 
many aspects traditionally treated by geophysical and space sciences. In 
this particular case, as we will see, it is especially true for the even zonal 
harmonic coefficients Jg of the terrestrial gravitational potential and their 
secular variations Ji. 

Recently, a test of the Lense-Thirring effect on the orbit of a test particle 
has been performed by Ciufolini and Pavlis (2004). They analyzed the data 
of the laser-ranged LAGEOS and LAGEOS II satellites in the gravitational 
field of the Earth by using an observable proposed in Iorio and Morea (2004). 
The data analysis spans 11 years; the 2nd generation GRACE-only EIGEN- 
GRACE02S Earth gravity model (Reigber et al. 2005a) released by the 
GeoForschungsZentrum (GFZ), Potsdam was adopted. The total accuracy 
claimed by Ciufolini and Pavlis is 5-10% at 1-3 sigma, respectively, but 
such estimate is controversial (Iorio 2005a) for various reasons. One of the 
major critical points is the impact of the systematic error induced by the 
secular variations Jt of the even zonal harmonics. Another critical remark 
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is that only one Earth gravity model was used by Ciufolini and Pavlis in 
their analysis. 

In this paper, we quantitatively investigate these issues. It will be shown 
that the error budget analysis performed by Ciufolini and Pavlis is opti- 
mistic. Our evaluation of the total error in the Lense-Thirring measurement 
ranges from 19% to 24% at 1-sigma level. 



2 The LAGEOS-LAGEOS II J 2 -free combination 

The adopted observable is the following combination of the residuals 5^ bs 
of the rates of the nodes of LAGEOS and LAGEOS II (Iorio and Morea 
2004) 

rrSLAGEOS i rrSLAGEOS II „ (a\ 

where c\ ~ 0.546, Slt = 48.1 mas yr _1 is the slope of the secular trend 
according to GRT, and //lt is equal to 1 in Einsteinian theory and in 
Newtonian mechanics. For the explicit expression of c\ see below Eq. (|7|. 
The combination of Eq. @ has been built up in order to cancel out the first 
even zonal harmonic J2 , along with its time- varying part which also includes 
its secular variations J2. The impact of the non-gravitational perturbations 
on the nodes of the LAGEOS satellites is of the order of ~ 1% of the Lense- 
Thirring effect (Lucchesi 2001; 2002; 2003; 2004; Lucchesi et al. 2004). 

The idea of using only the nodes of the LAGEOS satellites to disentangle 
the Lense-Thirring effect from J2 was proposed in Ries et al. (2003) in 
the context of the expected improvements in our knowledge of the Earth's 
gravity field from the dedicated GRACE mission. The linear combination 
approach was adopted for the first time by Ciufolini (1996) for his early, less 
precise tests with the nodes of the LAGEOS satellites and the perigee of 
LAGEOS II (Ciufolini et al. 1998). 

Let us write the observed residuals of the rate of the node of a satellite 
5Q in terms of the Lense-Thirring precession and of the classical precession 
induced by the Earth's quadrupole mass moment 

f^obs = ^lt/^lt + ^.2^2 + [higher degree terms], (5) 

where f^LT is as in Eq. Q and 

• d^geopotential 3 ( COsi f r\ 



in which R and n are the Earth's mean equatorial radius and the Keplerian 
mean motion n = \J GM/a 3 , respectively. The coefficients fl.e have been 
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explicitly worked out in Iorio (2003) up to degree I = 20; their numerical 
values for LAGEOS and LAGEOS II, whose orbital parameters are listed in 
Tabled can be found in Tabled If we write Eq. (JSJ) for both LAGEOS and 
LAGEOS II and solve for //lt the so-obtained linear system of two equations 
in the two unknowns /xlt and J2 it is possible to obtain Eq. © with 

C^AGEOS COS LAGEOS A" e LAGEOS I A ( OLAGEOSlA 

qLAGEOS II C os i L AGEOS II V 1 ~ e LAGEOS / V a LAGEOS / 

(7) 

Note that c\ is independent of the time span and of the even zonal harmonics: 
it is only fixed by the orbital parameters a, e, i of LAGEOS and LAGEOS 
II. 



3 The impact of the even zonal harmonics 

The combination of Eq. Q is, by construction, independent of the perturb- 
ing effects of degree t = 2 and order m = 0. In the case of J2 it can be 
checked by using the figures in Table El in order to calculate 

^LAGEOS + Ci ^LAGEOS II. (g) 

the result is zero. Note that the same also holds for the time-dependent 
perturbations like the I = 2, m = constituent of the 18.6-year tide. Indeed, 
by calculating the left-hand side of Eq. (j3J) with the values of the perturbing 
amplitudes of Table I and Table II in Iorio (2001b) it is possible to explicitly 
check this fact. 



3.1 The error due to the static part of the geopotential 

Eq. Q is affected by all the remaining even zonal harmonics of degree 
higher than two J4, Jq, Jg, along with their secular variations J4, Jq, .... 
The secular rate induced by J4, Jq, ... can be calculated as 

£ ^LAGEOS + Cl ^LAGEOS II) j £ (g) 

by using Table El for £l.e,£ = 4, 6, ... and the values for J^>4 of the chosen 
Earth gravity field model. Such aliasing shift is quite larger than the Lense- 
Thirring rate. Indeed, according to, e.g., the EIGEN-CG01C (Reigber et al. 
2005b) and EIGEN-CG03C (Forste et al. 2005) Earth gravity models, which 
combine data from CHAMP, GRACE and terrestrial measurements, the 



5 



nominal rate induced by the static part of the even zonals on the combination 
of Eq. © is ~ 10 5 mas yr^ 1 . It turns out that for the LAGEOS satellites 
the precessions due to the zonals with i > 16 can safely be neglected: indeed, 
their effect on the combination of Eq. (J1J amounts to 0.05 mas yr _1 ; the 
sensitivity is ~ 1 mas yr _1 . Thus, it is necessary to model the action of 
the low-degree even zonal part of the classical terrestrial gravitational field 
very accurately. To be more precise, in order to reduce the classical aliasing 
rates to a level ~ 1% of the Lense-Thirring effect it would be necessary to 
know J4 and Jq with an uncertainty of better than 2 x 10~ 12 and 4 x 10~ 12 , 
respectively. The present-day errors in J4 and Jq are, instead, 8.4 x 10 -12 
and 6.6 x 10~ 12 , according to the latest model EIGEN-CG03C. 

Unfortunately, it turns out that the present-day knowledge of the even 
zonal harmonics, represented by the latest models based on CHAMP and, 
especially, GRACE, is not yet good enough to allow for a reliable, model- 
independent 1% test of the Lense-Thirring effect with the combination of 
Eq. @. Indeed, at present GRACE seems to experience some difficulties in 
getting notable improvements in measuring the even zonal harmonics of low- 
degree (Wahr et al. 2004), contrary to the medium-high degree even zonals 
which, instead, do not pose problems to Eq. (J1J). As summarized in Table 01 
the systematic percent bias £//j^' potential induced by the mismodelled part of 
the geopotential on the gravitomagnetic effect ranges from ~ 4% for EIGEN- 
CG03C and EIGEN-GRACE02S to ~ 9% for GGM02S (Tapley et al. 2005). 
These figures are 1-sigma upper bounds obtained by linearly adding the 
absolute values of the individual mismodelled precessions, according to Eq. 
© with crj e>4 instead of J^>4, and then compared to the Lense-Thirring 
rate Slt- Since the signature of such source of bias is the same as that of 
the Lense-Thirring effect, it is not possible to fit and remove it from the 
time series without also affecting the relativistic signal of interest. It is 
only possible to evaluate as accurately as possible its aliasing impact on the 
Lense-Thirring trend. 

3.2 The error due to the secular variations of the even zonal 
harmonics 

Another source of non-negligible systematic error is represented by the sec- 
ular changes Jg>4 of the even zonal harmonics of low-degree (Cheng et al. 
1997; Bianco et al. 1998; Cox and Chao 2002; Dickey et al. 2002; Cox et al. 
2003; Cheng and Tapley, 2004). They are, at present, rather poorly known, 
as can be inferred from, e.g., Table 1 by Cox et al. (2003). Such temporal 
variations in the even zonals are mainly related to the Earth's lower mantle 
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viscosity features (Ivins et al. 1993). This topic has recently received atten- 
tion mainly due to the observed inversion of the rate of change of the Earth's 
quadrupole mass moment coefficient J2 which, since 1998, began increasing 
(Cox and Chao 2002). It is not yet clear if such an effect is a long-term 
feature or is short-term in nature; however, J2 is not relevant in the present 
context because Eq. Q is not sensitive to it. However, the possibility of 
interannual variations of J4 and Jq, depending also on the time span of the 
particular solution adopted, cannot be ruled out; they would affect, in prin- 
ciple, Eq. Q). In Tablc0]we quote the weighted means of the best estimates 
of Table 1 by Cox et al. (2003) for Ji>^. The uncertainties o" j f>4 reported 

in Table |H are the variances 1/c 2 = ^(l/of) of the distributions of the 
formal errors ai of Table 1 by Cox et al. (2003). Note that the so obtained 
a i are smaller than such formal errors. 

■Jl>4 

In the following we will perform quantitative evaluations of the impact 
of Ji>4 on Eq. Q. As a first, preliminary approximation, we will use the 
results of Table 0] which smooth out the possible interannual variations in 
the even zonals of interest. 

The effect of Jg integrated over the time grows quadratically. A first, a 
priori, quantitative evaluation of the aliasing impact of J4, Jq on the recovery 
of the Lense-Thirring effect can be performed by calculating the following 
quantity over a chosen observational time span T b s 

^ LA -^ • (10) 

For T bs = 11 years, it turns out that, according to Table El their impact 
amounts to ~ 12%. Note that since J2 does not affect the combination of 
Eq. Q and the errors in J4 and Jq are of the same order of magnitude of 
the nominal values, the situation does not change if we calculate Eq. (jTU|) 
with the a ; instead of their best estimates. 



4 Numerical simulations 

In this Section we investigate the impact of crj e and a j by means of numer- 
ical analyses. 



4.1 The simulated data 

Following Pavlis and Iorio (2002), we simulate with MATLAB a time-series 
of the combined node residuals of the LAGEOS satellites data (called Input 



7 



Model, IM) for the combination of Eq. @ by including in it 

• LT= Surt. Lense-Thirring trend as predicted by GRT according to 
Eq. @ 

. ZONDOT= Zti r t (^ AGEOS + Cl r2L AGEOS n ) (^) t 2 . Quadratic 

term due to the Ji according to Table |1J The numbers ri are ran- 
domly generated from a normal distribution with mean zero and unit 
standard deviation. 

• ZONALS= pSur-^t. Linear trend with a slope of x% of the Lense- 
Thirring according to Table El The number p is randomly generated 
as r£. 

• TIDE= a c OA c cos [(^r) t + / c ] Set of various tidal perturbations of 
known periods P. For the impact of such kind of perturbations on the 
orbits of the LAGEOS satellites see Iorio (2001b). The numbers a c , f c 
are randomly generated as p and rg. The tidal constituents included 
are Ki, Pi, S2, 165.565 and the 18.6-year tide 

• NOISE. We also included the effect of various other sources of non- 
systematic errors by using a generator of random numbers for all the 
points of the simulated time series. The resulting noise was a gaussian 
one and the amplitude could be varied. 

In a nutshell 

IM = LT + ZONDOT + ZONALS + TIDE + NOISE. (11) 

We include in our model the possibility of varying the length of the time 
series T b s , the temporal step At which simulates the orbital arc length, the 
amplitude of the noise and of the mismodelling in the perturbations and the 
initial phases of the sinusoidal terms in order to simulate different initial 
conditions and uncertainties in the dynamical force models of the orbital 
processors. The so-built IM represents the basis of our subsequent analyses. 

In order to check the reliability of the adopted procedure we first try 
to obtain some of the quantitative features of the real signal as released by 
Ciufolini and Pavlis (2004). For example, they fitted the 'raw' curve of Fig. 
2 (a) of their paper with a straight line only obtaining a Root Mean Square 
(RMS) of the post-fit residuals of 15 mas. To this aim we perform 5000 
runs for T obs = 11 years and x = 4 (EIGEN-GRACE02S) by fitting the 
simulated signals with a straight line (LF) clq + St. We calculate the RMS 
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of the residuals about the straight line fit. The averaged RMS amounts to 
15.6 mas. In Fig. ^ we plot, among other things, the simulated time series 
and the fitted straight line for one such run representing a generic set of 
initial conditions and mismodelling in the force models. It turns out that 
the RMS for such IM is 14.8 mas. This shows that our strategy represents 
a good starting point. 

The recovery of the slope S with LF is reliable because over 5000 runs 
the averaged error amounts to 

\ T obB =ll years 

nsu > 100) =0.8%. (12) 

5000 runs 



Slf / / 



4.2 Linear and quadratic fits 

In regard to the determined slope of the fitted data, a departure of ~ 7% 
from the predicted Lense-Thirring effect, which is present in IM, occurs for 
x = 4 (EIGEN-GRACE02S), i.e. 



Slf — Slt 



Slt 



\ T ohs =ll years 

100 ) ~7%. (13) 

/ 5000 runs 



The same holds for x = 3.8 (EIGEN-CG03C). For x = 6.3 (EIGEN-CG01C) 
and x = 8.7 (GGM02S) the departures amount to ~ 8% and ~ 10%, respec- 
tively. Instead, Ciufolini and Pavlis (2004) claim that the slope of their linear 
fit amounts to 47.4 mas yr _1 , i.e. a ~ 1% departure from the predictions of 
GRT. 

In the same set of 5000 runs, we also fit the basic IM, i.e. without any 
rescaling, with a quadratic polynomial (QF) ao + St + Qt 2 , finding 



5qf — Slt 



•Slt 



\ T b B =U years 

100 \ ~ 14%. (14) 

/ 5000 runs 



In order to assess the impact of the secular variations of the even zonal 
harmonics, we make the following preliminary remark. From a simple visual 
inspection of the plot of the simulated IM in Fig. ^ which refers to a IM 
built with l-o" j , it is not possible to infer anything about the presence or 
not of a quadratic term attributable to the secular variations of the even 
zonal harmonics. The parabolic signal is, indeed, very smooth and tends 
to affect the linear signal, especially as far as the early data are concerned. 
To be more quantitative, we assess the systematic error due to the rates 
of the even zonals by calculating for every run the difference between the 
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slopes obtained with the previously described LF and QF. Note that we do 
not include the sinusoidal perturbations in the fits because they are a minor 
problem. For example, the I = 2, m = 18.6-year tide does not affect Eq. 
(jlj) and over 11 years all the most relevant perturbations describe many full 
cycles; the longest period is that of the solar K\ tide perturbation on the 
node of LAGEOS, i.e. 2.85 years. Note that the approach of comparing 
the slopes of the fits performed with and without the features of which we 
want to assess the effect is the same adopted by Pavlis and Iorio (2002) and 
Ciufolini and Pavlis (2004) for the time-dependent harmonic perturbations. 
However, the problem of the secular rates of the even zonals was not treated 
in Ciufolini and Pavlis (2004), which only dealt with the time-dependent 
sinusoidal perturbations finding a 2% maximum error. 



Such results depend on the size of the parabolic signal introduced in the 
simulated IM. Indeed, if we conservatively triple the a j of Table 0] we get 
a discrepancy of 15%. The complete results of our simulations are shown in 
Tables El El and [7] for 1-Cj f , 3-<7j and 5-crj^, respectively. In Fig. [I]we plot 
the linear components of LF and QF for one of such runs. 

The same procedure has been repeated by rescaling the data of IM in 
order to shift the zero point of the time series in the middle of the data span. 
As a result, for T Q b s = 11 years the slopes of LF remain unchanged, showing 
the same 7 — 10% discrepancy with respect to the GTR value as before 
according to the different Earth gravity models used, while the slopes of QF 
do change in such a way that they are now identical to those of LF. In this 
case the magnitude of a j adopted in the IM does not influence the results. 
In Fig. |21 we plot the rescaled IM, the Lense-Thirring trend predicted by 
GTR and the linear components of LF and QF for one run. The results for 
different time spans and Earth's gravity models are shown in Table |H1 

5 Discussions and conclusion 

In this paper we addressed the problem of a quantitative evaluation of the 
impact of the current lingering uncertainty in the static and varying parts 
Ji and Ji of the even zonal harmonics of the terrestrial gravitational field on 
the measurement of the general relativistic Lense-Thirring effect with the 
combination of Eq. Q) involving the nodes of LAGEOS and LAGEOS II. 



We get 




(15) 
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To this aim, we realistically simulated the time series of the LAGEOS 
satellites data and we performed 5000 runs randomly varying the noise level, 
the initial conditions and the mismodelling of the included dynamical effects 
within the range of the currently known uncertainties. In each of such runs 
we fitted the basic simulated signal with a straight line and with a quadratic 
polynomial. The same procedure was also repeated by manipulating the 
simulated data rescaling them so that the zero point falls in the middle of 
the data set. 

5.1 The total error budget 

We propose the following error budgets for a data span of 11 years (all the 
obtained figures are 1-sigma bounds) 

• From a-priori analytical calculation (EIGEN-GRACE02S) 

Sfj^ al = ^(grav) 2 + (nongrav) 2 ~ y/(4 + 12 + 2) 2 + 2 2 = 18%. 

(16) 

The same also holds for EIGEN-CG03C. For EIGEN-CG01C we have 



V L ^ al = V(grav) 2 + (nongrav) 2 ~ ^{6.3 + 12 + 2) 2 + 4 2 = 20%. 

(17) 

GGM02S yields 

<5//$ al = V(grav) 2 + (nongrav) 2 ~ ^(8.7 + 12 + 2) 2 + 2 2 = 23%. 

(18) 

We have included in our evaluations also the errors due to the tides 
(2%) and the non-gravitational perturbations (2%). 

• From numerical simulations and fits of the basic simulated data (EIGEN- 
GRACE02S) 

Suffi 1 = V(grav) 2 + (nongrav) 2 ~ ^/(7 + 13 + 2) 2 + 2 2 = 22%. 

(19) 

The same also holds for EIGEN-CG03C. For EIGEN-CG01C we have 



Sfiffi 1 = V(grav) 2 + (nongrav) 2 ~ y/(8 + 13 + 2) 2 + 2 2 = 23%. 

(20) 

GGM02S yields 

(SMut" 1 = V(grav) 2 + (nongrav) 2 ~ ^(10 + 13 + 2) 2 + 2 2 = 25%. 

(21) 
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• From numerical simulations and fits of the rescaled simulated data 
(EIGEN-GRACE02S) 

<5/x£f 1 = V(grav) 2 + (nongrav) 2 ~ x /(7 + 2) 2 + 2 2 = 9%. (22) 

The same also holds for EIGEN-CG03C. For EIGEN-CG01C we have 

6^ = ^/(grav) 2 + (nongrav) 2 ~ ^(8 + 2) 2 + 2 2 = 10%. (23) 
GGM02S yields 

6^ = ^(grav) 2 + (nongrav) 2 ~ y/(W + 2) 2 + 2 2 = 12%. (24) 

In regard to the evaluation of the error of gravitational origin, the secular 
variations of the even zonal harmonics were not solved for in the latest 
Earth gravity field models based on CHAMP and GRACE: they were held 
fixed to default values obtained from multi-year data analysis of the data 
of the constellation of laser-ranged geodetic satellites. Thus, the recovered 
values for Ji are not independent of J%. Thus, caution advices to linearly 
sum the biases of the static part of the geopotential and of their secular 
variations in assessing the total systematic error of gravitational origin in 
the Lense-Thirring test. 

Instead, Ciufolini and Pavlis (2004) at the end of the Section Total un- 
certainty, pag. 960, and in the Supplementary Information .doc file add 
in quadrature the doubled error due to the static part of the geopotential 
(i.e. 2x4% value obtained from the sum of the individual error terms), their 
optimistic evaluation of the error due to the time dependent part of the 
Earth gravity field (2%, not doubled) and the non-gravitational error (2%, 
not doubled) getting V8 2 + 4 2 + 4 2 % = 10%. On the other hand, in the 
Supplementary Information .doc file it seems that they triple the 3% error 
due to the static part of the even zonal harmonics obtained optimistically 
with a root-sum-square calculation and add it in quadrature to the other 
(not tripled) errors getting V9 2 + 2 2 + 2 2 % < 10%. 

These considerations and our results show that the 1-sigma 5% total 
error claimed by Ciufolini and Pavlis (2004), who used only one Earth's 
gravity model, is optimistic. Moreover, the present-day uncertainty in the 
knowledge of the static part of the geopotential does not yet allow for a 
fully reliable, model-independent determination of the gravitomagnetic ef- 
fect with the combination of Eq. @. It should also be considered that the 
evaluations presented here might turn out to be optimistic because they are 
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based on the smoothed values of Table |1J The impact of possible interannual 
variations of J4 and Jq, depending also on the particular observational time 
span considered, has not been considered here and should deserve further 
investigations. 

An inspection of Tables 00 shows that the secular variations of the even 
zonal harmonics may represent a limiting factor also for longer and shorter 
time spans. Unfortunately, the expected improvements in the knowledge of 
the static part of the geopotential from the forthcoming GRACE solutions 
will not ameliorate the situation with respect to the problem of the j£. 

5.2 How to improve the accuracy and the reliability of the 
Lense-Thirring tests 

Among possible alternatives the most promising one is the launch of an- 
other satellite of LAGEOS-type, like LARES or OPTIS (Lammerzahl et al. 
2004). As shown in Iorio (2005b), the improvements in our knowledge of 
the gravitational field from the GRACE mission would allow to make the 
requirements on the orbital geometry of such a new laser target much less 
stringent than in the originally proposed LARES mission. In particular, by 
combining the data of the new spacecraft with those of LAGEOS and LA- 
GEOS II, it would be possible to reduce its semimajor axis from 12270 km 
to 7500-8000 km: this would allow to greatly reduce the costs of the mis- 
sion. Also the inclination could experience large departures from the value 
of LARES (i = 70 deg) without affecting the outcome of the experiment. 
The accuracy in measuring the Lense-Thirring effect with a three nodes 
combination including also LAGEOS and LAGEOS II is ~ 1% at 1-sigma 
level. 
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Table 1: Orbital parameters of the existing LAGEOS and LAGEOS II satel- 
lites and their Lense-Thirring node precessions in mas yr _1 . 



Satellite 


a (km) 


e 


i (deg) 


£1 L t (mas y 




LAGEOS 


12270 


0.0045 


110 


31 




LAGEOS II 


12163 


0.0135 


52.64 


31.5 
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Table 2: Coefficients of the node classical even zonal precessions of 
LAGEOS and LAGEOS II up to degree t = 20 in mas yr" 1 . 



1 


LAGEOS 


LAGEOS II 


2 


4.191586788514 x 10 11 


-7.669274920758 x lO 11 


4 


1.544030247472 x 10 11 


-5.58637864293 x 10 10 


6 


3.25092246054 x 10 10 


4.99185703735 x 10 10 


8 


2.1343038821 x 10 9 


1.10707933989 x 10 10 


10 


-1.4885315218 x 10 9 


-2.2176133068 x 10 9 


12 


-7.703165634 x 10 8 


-1.1555006405 x 10 9 


14 


-2.097322521 x 10 8 


2.5803602 x 10 6 


16 


-3.04891722 x 10 7 


8.81906969 x 10 7 


18 


2.7037212 x 10 6 


1.25437446 x 10 7 


20 


3.3458376 x 10 6 


-4.8988704 x 10 6 



Table 3: Uncertainty of the systematic error 5fif^ po due to the static part 
of the even zonal harmonics on the Lense-Thirring effect for the combination 
of Eq. according to the variance matrices of the Earth gravity models 
EIGEN-CG03C, EIGEN-GRACE02S , EIGEN-CG01C and GGM02S. The 
linear sum of the absolute values of the individual terms, according to Eq. 
© with aj e instead of Ji, have been used in order to give realistic 1-sigma 
upper bounds. 



EIGEN-CG03C 


EIGEN-GRACE02S 


EIGEN-CG01C 


GGM02S 


fyg? pot 3.8% 


4.1% 


6.3% 


8.7% 
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Table 4: The reported values for J2, J4 and Jq are the weighted means of the 
best estimates of Table 1 by Cox et al. (2003) in units of 10~ n yr" 1 . The 
uncertainties a j are the variances of the distributions of the formal errors 
of Table 1 by Cox et al. (2003) in the same units. 



I = 2 


£ = 4 


e = 6 


J e -2.113 


-0.6992 


-0.3594 


at 0.0810 
il 


0.2029 


0.1765 



Table 5: Results of our tests for the evaluation of the 1-sigma level impact 
of Je on the measurement of the Lense-Thirring effect with the combination 
of Eq. (J3J. No rescaling has been applied to the simulated time series. 
Alf-qf is the averaged difference between the slopes of the fitted linear 
trends in the linear fit-only (LF) and the quadratic fit (QF). Alf-lt is the 
averaged difference between the slopes of the fitted linear trend in the linear 
fit-only (LF) and the Lense-Thirring trend as predicted by GRT. Aq F _lt is 
the averaged difference between the slopes of the linear part of the quadratic 
fit (QF) and the Lense-Thirring trend as predicted by GRT. A time step of 
At = 15 days and EIGEN-GRACE02S have been used. 



T obs (yr) Alf-qf (%) A LF _ LT (%) A QF _ LT (%) 



3 


22 


11 


24 


6 


10 


8 


14 


11 


13 


7 


14 


20 


12 


7 


11 


30 


9 


9 


4 



Table 6: Results of our tests for the evaluation of the 3-sigma level impact 
of Je on the measurement of the Lense-Thirring effect with the not-rescaled 
simulated time series. A time step of At = 15 days and EIGEN-GRACE02S 
have been used. 



T obs (yr) Alf-qf (%) A LF _ LT (%) A QF _ LT (%) 



3 


22 


10 


25 


6 


11 


9 


14 


11 


16 


10 


14 


20 


18 


14 


11 


30 


20 


20 


4 
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Table 7: Results of our tests for the evaluation of the 5-sigma level impact 
of Jn on the measurement of the Lense-Thirring effect with the not-rescaled 
simulated time series. A time step of At = 15 days and EIGEN-GRACE02S 
have been used. 



T obs (yr) Alf-qf (%) A LF _ LT (%) A QF _ LT (%) 



3 


22 


11 


24 


6 


12 


11 


14 


11 


18 


14 


14 


20 


25 


23 


11 


30 


33 


33 


5 



Table 8: Results of our tests with the rescaled simulated time series. Alf-qf 
is now zero and Alf-lt and Aqf-lt are equal and denoted as A (in per- 
cent). There is no dependence on the size of a j . A time step of At = 15 
days and different Earth's gravity models have been used. 



T obs (yr) 


A (EIGEN-GRACE02S) 


A (EIGEN-CG01C) 


A (GGM02S) 


3 


10 


12 


13 


6 


9 


10 


11 


11 


7 


8 


10 


20 


4 


5 


7 


30 


3 


5 


7 
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Figure 1: Basic, not-rescaled simulated time series (•), predicted Lense- 
Thirring trend (— ), trend of the straight line-only fit (— — ) and trend of 
the quadratic fit (• — •) for T Q b s = 11 years, At = 15 days and EIGEN- 
GRACE02S. The effect of Ji is present at 1-sigma level, according to Table 
0J in the simulated time series. The difference between the slopes of the 
two fitted linear trends amount to 13% of the Lense-Thirring effect. A 7% 
discrepancy between the straight-line only fit and the predicted LT slope 
occurs. See Table 03 
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Figure 2: Rescaled simulated time series (•), predicted Lense-Thirring trend 
(— ), trend of the straight line-only fit (— — ) and trend of the quadratic fit 
(• - •) for T obs = 11 years, At = 15 days and EIGEN-GRACE02S. The 3 t 
effect is present at 1-sigma level, according to Table 01 in the simulated time 
series. The difference between the slopes of the two fitted linear trends is 
now negligible. A 7% discrepancy with the predicted LT slope occurs. See 
Table El 
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